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In this paper, analytical formulas for elastic neutrino-fermion scattering cross sections are presented. 
These are evaluated in the limit where the neutrino mass is neglible compared to the incoming 
momentum. The fermion mass is taken into account, however, yielding a good approximation to 
realistic scattering reactions. The intended level of readership is for undergraduates/graduates in 
the fields of relativistic quantum mechanics and quantum field theory. In addition to a detailed 
and pedagogical derivation of the analytical expressions, numerical results are also provided for 
/ = e,n,p. 

PACS numbers: 13.15.+K, 25.30.Pt 



I. INTRODUCTION 

The motivation for writing a paper on analytical formulas for elastic neutrino-fermion scattering cross sections is 
that i) their explicit derivation is hard to find in literature in general, ii) their derivation provides very useful insight 
into how particle reactions are analyzed in the framework of quantum field theory and Feynman diagram analysis, 
and Hi) that some important equations and aspects of the analysis are often left out or sparsely treated in original 
textbooks. This paper should be of particular interest to students looking for a detailed review of how scattering 
reactions are dealt with, in addition to some important observations and techniques involved in this task. Textbooks 
such as Mandl&Shaw [l[ provide excellent introductory texts to the topic of quantum field theory, and equips the 
reader with many tools to handle scattering reactions. The "dirty" details of such calculations can be hard to find, 
however, and are highly valuable and appreciated by students trying to understand the bits and pieces of how one 
arrives at a given result. 

This paper is organized as follows. In Sec. [n] the concept of a reaction cross section is briefly reviewed to set the 
notation. Sec. HHI takes on neutrino interactions in electro-weak theory, which is the necessary framework for studying 
neutrino-fermion scattering. To set the notation, the relevant Feynman symbols are presented in Sec. IIVI The general 
expressions for neutrino-fermion cross sections are derived in Sec. [V] and a summary is given in Sec. IVII Unless 
specifically stated otherwise, natural units h = c = 1 will be used. 

II. CROSS SECTION 

The differential cross section da says something about how incoming particles are scattered in space, and the 
likelihood of the scattering to occur at all. More precisely, since the interpretation of a is the number of particles 
scattered into a specific final state per unit time divided on the incident flux, the expression for da is 



wV yj Vdp' f 

where w is the transition probability from \i) to |/) per unit time, v T0 \ is the relative velocity between incident and 
target particles, and Y[f (2^3 is the number of states with momentum in the interval (p^,py + dp^-). In the special 
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case of two incoming particles, the differential cross section reads 



(2) 



where pi are the 4-momenta of the incoming particles, p'j are the 4-momenta of the outgoing particles, and M. is the 
Feynman amplitude of the process. More on this in Sec. IHII 



Calculations are often simplified when considering the process from the Center-of-Mass (CM) system. Here, Eq. 
is reduced to 



da 



iP'il 



^dfi/cM 64tt 2 (E 1 + E 2 ) 2 |pj 
For a detailed derivation of these quantities, consider for instance Ref. [|. 



(3) 



III. NEUTRINO INTERACTIONS IN ELECTRO- WEAK THEORY 

In the following, we shall not consider neutrino interaction with the electromagnetic field A^ix), since the electromag- 
netic transition moments of the neutrino are too weak for any significant coupling to photons. The upper limit for 
the neutrino magnetic moment is ~ 10~ 10 /Lt B , taken from the Particle Data Group |2|. When dealing with neutrino 
interactions it is sufficient to consider the Lagrangian densities involving W and Z exchange for reactions involving 
V\ . We shall stick with the notation of Ref. [l| . Now, the interaction term reads 
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(x)Wt(x) + J^(x)W fi (x)] 
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cos 8w 



J%(x)Z°(x), 



(4) 



where j!L{x) and are the charged and neutral currents 



J^(x) = ^l(x)^(i- 75 )n(x) 
1 



(5) 



In Eq. JSJ, l(x) and vi{x) represent the second quantized lepton fields. Furthermore, g % v and g\ are weak coupling 
constants, 9w is the Weinberg angle, and 7^ are the Dirac gamma matrices. Wp(x) and Z®(x) describe the vector 
boson propagators found in Fig. [5] Eq. J2J| is valid for all neutrino energies, but we would like to derive an expression 
for the low-energy limit which often occurs in realistic scattering experiments. More specifically, consider the limit of 
the external momenta being much smaller than the vector boson masses mw and mz- This means that 



lim 



lim 

m 2 »fc 2 



m 



w 



k 2 — m 2 ? + ie 
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(6) 



with g^u = diag(l,-l,-l,-l). Making these substitutions corresponds to the interaction transformation in Fig. ^ and 
thus provides us with the effective interaction Lagrangian density for neutrino interactions 



r cff_ r eff , r eS. _ ( 9 \ 2 1 T n r \ ( 



cos 9w 



— J M J f 



FIG. 1: Effect of fow-energy approximation with dominant mw and mz- 



The gauge boson masses are related through the Weinberg angle 9w by mw / mz = | cos 9w \ ■ u and d quarks can also 
be incorporated into this model, merely by adding extra terms to the currents in Eq. JSJ , thus producing 

J w = - 75)^+^7^(1 - 7sH 

i 

J z = \ E [^W - fi-ft)** + l-ftfv - 9aJ 5 )1 
I 

+ ur(9v ~ 9a75)u + drf<j& - g d A ^)d] , (7) 

where dg = dcosOc + ssin9c- The Cabbibo angle 9c has been experimentally determined to cos 9c ~ 0.98, so we 
shall exclude the s quark part from now on and set dg = d. The introduction of quark terms enable us to describe 
processes such as v e + n — > p + e. When adding the quark terms, we have used the general expression for the neutral 
current found in Ref. namely 



J z = \ E ^ fe 1 " 75) - 2Q* sin 2 9 W 



(8) 



where i = (I, vi, u, d), while If is the belonging particle isospin and Qi is the particle charge in units of e. 



Now, the Hamiltonian density TL is related to its Lagrangian density £ by TL — jnfcv — £, with <fi = <j)(x) denoting 
the various fields appearing in £, i.e. l(x),vi(x) et.c. The second quantized expansion for Dirac fields reads 



^) = E 



2VE n 



~K*) = E 
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a s (p)u s (p)e-^ + &t(p) Us (p) e ^ 
6 s (p)U s (p)e- i ^+ a t(p) IIs (p)e^ 



(9) 



and similar for Ui(x),u(x),d(x). Here, E p is the energy of the lepton with momentum p, {a s (p), a|(p)} and 
{b s (p), i>J(p)} are annihilation and creation operators for particles and antiparticles, respectively, while {w s (p), w s (p)} 
are Dirac spinors. Note that l(x) = ^(x)^ . 

Constructing the interaction Hamiltonian density for a process is a matter of straight-forward calculation, which is 
considerably simplified by making the observation that C does not contain any time derivatives of the fields, i.e. 
C = -U. We arrive at H cS = Hf s w + Hf z , where Hf w = % Jw»J$, Uf z = ^J%J Zt i- The Fermi constant is 
defined as Gp — V2g 2 /8m^,. 



IV. FEYNMAN RULES 



The Feynman diagram of a process is an indispensible tool for analysis of any particle reaction, in order to distuingish 
between different outcomes and to equip us with an elegant manner of calculating important physical quantities, such 
as the cross section. The Feynman amplitude can be read from the diagram with some training and quantifies the 
process. For our needs, we shall only give a short overview of the symbols and algebraic equivalents used in Feynman 
diagrams describing weak processes to set the notation. This is seen in Fig. |2 
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EXTERNAL LINES SYMBOL FACTOR IN M 

(P,r) 

Initial lepton line / or vi ^ ^ u r (jp) 

Initial lepton line l + or vi ^ ^ w r (p) 

Final lepton line l~ or f ( ' u r (p) 

• ► • 

(p, r) 

Final lepton line l + or z/; v r (p) 



PROPAGATORS SYMBOL FACTOR IN M 

k 

Internal W boson ^\f\J\f\f\J\j K-9°g+fc° W™w) 

k 

Internal Z boson * \J \J \J \J \J \j k'-m' z +ie 



CHARGED CURRENT VERTEX SYMBOL FACTOR IN M 



n a 

Interaction between leptons {I, v{\ ». ► ^v§^ Q ^ — 

and W, I G {e, /Li, t} 

W 



NEUTRAL CURRENT VERTEX SYMBOL FACTOR IN M 



i a i 



Interaction between particle i and Z, 
i e {l,vi,u,d},l G {e,^,r} 




l a (9v -9a15) 



FIG. 2: Algebraic expressions for Feynman symbols in standard weak theory, relevant for our needs. 



V. NEUTRINO-FERMION SCATTERING 

In this section, we first give attention to neutral current neutrino-fermion scattering, i.e. mediated by the Z° boson. 
Then, we generalize to reactions mediated by the W boson. Then it is time to say something about the coupling 
constants appearing in the expressions for the total cross section. 
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A. Neutral current scattering 

To begin with, we shall consider Vij (I ^ /) scattering. This process is mediated only by the Z° boson, if we disregard 
the Higgs boson as a gauge particle. This is justified by the fact that the Feynman amplitude corresponding to the 
scattering process mediated by H is of 0(m Vl m t / 'm 2 H ) compared to the Z° diagram. The only relevant lowest order 
Feynman diagram for vif (I ^ f) scattering is shown in Fig. [3] 

n(pi) Mp'i) 



Z° y kz = pi - p[ 



/(Pa) /«) 
FIG. 3: Lowest contributing order Feynman diagrams for Vi + / — > vi + /. 

The Feynman amplitude is given by 

M = U - {1 )h 2c^ 7 (W-ffi75)K(l)i( fc| _ m| + i£ ) 

Observe that we have shortened down u Si (pi) to u(i). 

These expressions can be simplified a great deal, when making some considerations. It is reasonable to expect 
mly ^> k^y when using mw = 80.4 GeV for the intermediate vector boson W^. The same inequality goes for 
mz = 91.2 GeV and mu > 112 GeV. By applying this low-energy limit, we regain the effective propagator from Sec. 
IIIII arriving at 

M = ^|£u v ,(l')[ 7 Q (^ -g^ b )]u Vl (l)u0)[ la {g v -g f Al ,)]u { {2). (10) 

We now pursue the absolute square of the Feynman amplitude, |.M| 2 , which is part of the differential cross section 
from Sec. [HI In a perfect world, the polarization, i.e. spins, of the reacting particles would be known. Unfortunately, 
this is the exception rather than the rule. If one does not know the initial and final spins of the neutrinos and 
electrons, this must be compensated for by averaging over the incoming spin states and adding the final states. This 
normally corresponds to an averaging factor 1/2 for each incoming particle. 



Neutrinos, however, are special cases. Although they are massive, their mass is so small that it is ordinarily set to 
zero in the standard model. This means that there is only one helicity state for the neutrino, since it consequently 
travels with the speed of light and never can be overcome in any Lorentz frame. Thus, there is only one spin state to 
consider for the neutrinos which leads to the substitution 

\m\^\Y,\m? 

Si,s' f 

in order to obtain the unpolarized cross-section. By exploiting the fact that lA^I 2 = M.M.* and {^y — the 
problem is attacked. Inserting Eq. (|10|l gives 



^^(l') 7 Q (ffy ! " ff^75)^ / (l)w/(2 , )7 Q (5^ " ff£ys)«/(2) 

-^■"/ ( 2 )7/8 (5v -5A75)w/(2')^(l)7 /3 (5y -3^75)^(1') 

We stride on by recalling the completeness relations for the Dirac spinors it and w 

y^M; a (p,g)^(p } g) = {$ + mi) a p, ^2vi a (p,s)vip(p,s) = {$-mi) a p. 

s s 

Attaching indices on spinors and matrices provides us with 
1 



£|M| 2 = G|]T 



We designate 



Uu lK (l') -fitf-sKis) u UlX (l)u fj {2') i a (g f v - g f A i 5 ) u„(2) 

L J yd 



K A 



X «/e(2) 7/3(5y -3a75) u/t (2>^,t(1) 7 /3 (ffy " fl^Tto) u VlP {\') 

L J er L J op 

xTr | (#2 + m f)l a {g f v - 75)(i$2 + m f ) lf] {g f v - 5^7s)| ■ 

® = Tr{ + m 1/( ) 7 «(^ - 3^5) (A + m Vl )^{g$ - <$7s) \ 



® = Tr j(#> + m fh a (g{, -3a75)(^2 + m f )^ {g 1 v ~g A To)j, 

and consider each of the traces separately. Using the general relations for traces 
Tr{A + B + ...} = Tr{A} + Tr{B} + ... 
Tr{ 7 5 } = Tr{ 7 5 7 Q } = Tr{ 7 5 7 Q 7 /3 } = Trfr^VV} = 
TtfrVVVY} = _ 4ie a/3 7 <5 

Tr{ 7 a 7 ^... 7 ' 3 7 cr } = if 7 a 7 /3 ... 7 p 7 CT is an odd number of 7 -matrices, 
where £ a Pi s is a completely antisymmetric tensor, we obtain 

® =8pi^ \^g A l ^ avp + \wf + Otf ) a ]fa /i V /3 + <V a - rTv ap ) 



J J 



^gv9A e aapil + o 



in the limit m^, — > 0. From Eq. we see that the differential cross-section in the CM system must be 



X 



IPil 



where 



Vdfi/cM 647r 2 (^i + £ 2 ) 2 |j>i| 
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The product of the two traces is 

© x ® = 128<^ [Wi){pm) ~ (p'iP2)(pip' 2 )] 

+ 32[(g f v ) 2 + (^) 2 ][(.^) 2 + ( 9 m MP2)(P1P' 2 ) + (p^)(piP2) 
- 32m 2 [( 5 £) 2 - + Gtf ) a ](pipi)]. 

The situation now looks like Fig. ^in the CM frame with 4- momenta 

p x = (E,0,0,E), 
P2 = (yjE* + m 2 f ,0,0,-E), 
p[ = E(l, cos (j)sm0, sine/) sin 9, cos 9), 



p' 2 = (y E 2 + to 2 , — J5 cos sin 6*, — E'sin <^>sin0, — J5 cos ( 
where _E = E CM is the CM kinetic energy of the neutrino. 



(16) 



(17) 



FIG. 4: vi f (I f) scattering in the CM frame. 
Since all 4-momenta now are known, insertion into Eq. (|16J) gives 



© x ® = 128g v g A g^g A ' [(EjE 2 



e 2 ) 2 - (eje 2 



+ 32{(g v ) 2 + {g f A f][{g^f + (g^) 2 } [(E^jE 2 + m 2 +E 
- 32[( 3 £) 2 - (g A ) 2 }[(g^) 2 + {g'2) 2 ]E 2 (l - cos0)] . 



2\2 



E 2 cos6>) 2 l 



(EJE 2 



The total cross section for vif (I ^ /) scattering is now within our grasp, and reads 



E 2 cos6>) 2 



(18) 



*Wif (i /)] 



G% 



32n{E + JE 2 + to 2 ,) 2 Jo 



© x (2) sin 9d6. 



From Eq. (|18|) . it is clear that we must evaluate three types of integrals, namely 



sin0d0 = 2, / cos6»sin(9 AO = 0, / cos 2 9 sin 9 A9 = - . 

Jo Jo 3 



(19) 



(20) 



Note that s = (pi +P2) 2 = [E+ \J E 2 + m 2 ) 2 . At this point, we make a little "cheat" and reveal the neutrino coupling 

constants gy and g A to be 1/2. We shall remedy this in Sec. IV CI Inserting our integrals Eq. (|20|l reduces Eq. i|19|) 
to 



<y[vif (f + I)] = 



G 2 F (s - m 2 f 
Airs 



{9v + 9a? + {9v-9 A ) 



f _ J\2 



I 2 (s — TO 2 ) 5 



3.S 2 



U) 2 -( 9 f A ) 2 }^ 



(21) 
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The cross section for antineutrino-fermion scattering, i.e. T>if{f 7^ I), is obtained from the following argumentation. 
In the limit m Vl — > 0, 1 = e, /1, r, neutrinos are always left-handed while antineutrinos will be right-handed. Thus, 
only a parity transformation V on the neutrino vertex part of the Feynman amplitude Eq. IjlOl) is required to 
obtain 77;/' (/ ^ I') scattering. Now, a scalar product of two axial vectors is invariant under a parity transformation, 
and so is the scalar product of the vector quantities as well. Recall that axial vectors j/ M = (y°,y) transform as 



y^ 1 = (—y°,y) under parity, while vector quantities z M = (z°,z) transform as z h 



= (z°,-z). Thus, it 



follows that a scalar product of one axial and vector quantity is not invariant under parity transformations, since 



2/% 



y°z + yz 



(22) 



As a consequence, all mixed terms gyg\ l wm change sign. This means that the cross section for antineutrino-fermion 
scattering is found by making the substitutions 



(23) 



For vif (I 7^ /) scattering, g v g A wu l a l so change sign, yielding 



vplf (I * /)] = "Wlf (I * /)]■ 



(24) 



B. Charged current scattering 



Having derived the general expression for elastic neutrino-fermion scattering processes mediated by the Z° boson, we 
can now include the case / = /. This process is mediated by the W boson in addition to Z , thus providing Feynman 
diagrams as shown in Fig. [S] 



n(j>i) 



Z° f k z =pi- pi 



W f = Pi - P2 



i( P2 ) i{ P ' 2 ) *(p0 u '(p'i) 

FIG. 5: Lowest contributing order Feynman diagrams for i>i + I — » v\ + 1. 

The total Feynman amplitude in the low-energy limit k 2 <C {m 2 z , m^} is now A4 = A4z + M-w-, where 



Mi 



' h ° F -u Vl (!') [ 7 Q (<# - 7s)K (2') [ 7 a (g l v - (2) , 



V2 
-iG 



J^ff = -7=-^ (2') [l a (1 - 7&)K (IK (1') [ 7 a (1 - 1b)]ui (2). 
By using the Fierz identity 

"l[7 Q (! - 75)]ti2W3[7a(l - 75)]"4 = Ul[7 Q (l - 75)]w4U 3 [7a(l ~ 75)]«2 



(25) 



(26) 



for arbitrary spinors Ui, Mw obtains the same sequence of spinors as Mz- In addition, we take advantadge of the 
fact that gy = $^4 = 1/2, I — e, fi, t. This leaves us with 



7 Q (<# ~ ^'75)1^,(1)^ (2') k[G?y + 1) - (g l A + l)71s]l«i(2) 



(27) 
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This is exactly the amplitude Eq. (|10|) with the replacements g v — ► ijy + 1, g A — > <A + 1. Therefore, we are 
immediately able to write down the total cross section for elastic v\l scattering, namely 



aM = lim o[ Vl f (I ? /)]. 



(28) 



4-3^ + 1 



C. Coupling constants 

The derivation of the coupling constants starts from the condition of U(l) and SU(2) gauge invariance for the 
standard electro-weak Lagrangian. Although an interesting exercise in itself, it is outside the scope of this paper, 
and we simply state the results in Tab. [I] 



TABLE I: Relevant weak coupling constants for our needs with I = {e, n, r} from Ref. Q. 



Particle 


Axial coupling qa 


Vector coupling gv 


I 


-1/2 


2 sin 2 6w -1/2 


VI 


1/2 


1/2 


u 


1/2 


1/2 -4 sin 2 6 w /3 


d 


-1/2 


2 sin 2 6W/3 - 1/2 



It should be mentioned that our factor g A = 2g A + g d A is not exactly equal to 1/2. This is due to the axial current 
anomaly, which states that the conservation law for the axial current is broken. Ref. [4| gives a proper treatment of 
this topic. Neutron beta-decay experiments have determined the value of axial vector coupling constant for protons 
to g p A ~ 0.635 (see Ref. 0). 



VI. SUMMARY 



We have found that the elastic cross section for v\f scattering, where I = e,fx,T and / ^= I is a fcrmion with vector 
coupling strength g v and axial coupling strength g A , is given by 



°\nS (/ + 0] = 



G 2 F {s-m)f 



(9v+9 A ) 2 + (9v-9 J A ) 



m 2 f ) 2 



3s 2 



H) 2 - ( 9 f A ) 2 }^ 



(29) 



where s — (E CM + y E 2 M + m 2 ) 2 = 2E LAB nif + to 2 . For charged current v\l scattering, the cross section is obtained 
through 



a[ Vl l]= lim a[ H f(f¥=l)l 



(30) 



Two specific limits of Eq. 129(1 are of particular interest. By taking m/ — > 0, the resulting cross section corresponds 
to elastic neutrino scattering on light leptons. This approximation is very good for / = {e,/i}. In this case, Eq. 129|) 
reduces to 



lim ah/ (fjkl)] = 9bd 

m f —>0 07T 



9v+9a + 9v9a 



(31) 
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where we have dropped the supercript / for the fermions. Another interesting scenario is the limit £ LAB /m/ <C 1. 
The physical interpretation of such a limit is neutrino scattering on heavy fermions, e.g. protons, where the resulting 
recoil energy of the fermion is small. One finds that 



15m a[ H f(f^l)} = 

-BLAB/m/«l 



r>2 ZP2 
F LAB 



9 2 v 



(32) 



This result agrees with theoretical predictions and experimental data (see e.g. Ref. 0). Fig. [(^contains a plot of Eq. 
(|21|l for the most common fermion scattering components on Earth, / = e, n,p. It is also instructive to consider the 
mean free path length / m f p for the neutrinos in a medium. This quantity is given as i m f p = \/{NfOf), where Nf 
is the number density of scattering component /. Taking / = e, n,p, the resulting mean free path length is shown in 
Fig. [7| for varying mass densities p = ^2f Nfirif, assuming an electrically neutral medium with Nf = N . 



1.4 - 
1.2 - 



— 1 " 
E 




E [MeV] 



FIG. 6: The elastic neutrino-fermion scattering cross section Eq. (12 U for / = e,n,p as a function of the CM neutrino kinetic 
energy E. 
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FIG. 7: Log-plot of mean free path length Z mfp in meters as a function of CM neutrino kinetic energy E for various mass 
densities. 



